o 
o 

(N 

> : BARYOGENESIS 
O 

CSl ' James M. Cline 

(N ■ 



> 
in 

a^ 
o 
\o 
o 



X 



McGill University. Dept. ofPhy.sics. Montreal. Qc H3A 2T8. Canada 



1 



Photo: width 7.5cm height 11cm 



Contents 



1 . Observational evidence for tlie BAU 

2. Sakiiarov's Conditions for Baryogenesis 

2.1. B violation 

2.2. Loss of thermal equilibrium 

2.3. C, CP violation 

2.4. Examples of CP violation 

2.5. More history 

3. Example: GUT baryogenesis 
3.1. Washout Processes 

4. B and CP violation in the Standard Model 
4.1. CP violation in the SM 

5. Electroweak Phase Transition and Electroweak Baryogenesis 

5.1. Strength of the phase transition 

5.2. EWPT in the MSSM 

6. A model of Electroweak Baryogenesis: 2HDM 
6.1. EDM constraints 

7. EWBG in the MSSM 

8. Other mechanisms; Leptogenesis 
References 



3 



1. Observational evidence for the BAU 



In everyday life it is obvious that there is more matter than antimatter. In nature 
we see antimatter mainly in cosmic rays, e.g., p, whose flux expressed as a ratio 
to that of protons, is 

^ ~ 10"^ (1.1) 

P 

consistent with no ambient p's, only p's produced through high energy collisions 
with ordinary matter. On earth, we need to work very hard to produce and keep 
e"*" (as in the former LEP experiment) or p (as at the Tevatron). 

We can characterize the asymmetry between matter and antimatter in terms of 
the baryon-to-photon ratio 

rj^I^E^ (1.2) 
where 

Hb = number density of baryons 
ng = number density of antibaryons 
n-y = number density of photons 
C(3) , 

= — ^QifT , £f* = 2 spm polarizations (1.3) 

and C(3) = 1.20206 . . . This is a useful measure because it remains constant 
with the expansion of the universe, at least at late times. But at early times, and 
high temperatures, many heavy particles were in thermal equilibrium, which later 
annihilated to produce more photons but not baryons. In this case, the entropy 
density s is a better quantity to compare the baryon density to, and it is convenient 
to consider 

s 7.04 ' ^ ' 

where the conversion factor 1/7 is valid in the present universe, since the epoch 
when neutrinos went out of equilibrium and positrons annihilated. 
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Historically, 77 was determined using big bang nucleosynthesis. Abundances 
of ^He, ''He, D, ^Li and ^Li are sensitive to the value of 77; see the Particle Data 
Group [1] for a review. The theoretical predictions and experimental measure- 
ments are sunmiarized in figure 1, similar to a figure in [2]. The boxes represent 
the regions which are consistent with experimental determinations, showing that 
different measurements obtain somewhat different values. The smallest error bars 
are for the deuterium (D/H) abundance, giving 



for the two experimental determinations which are used. These values are con- 
sistent with the "'He abundance, though marginally inconsistent with ""Li (the 
"Lithium problem"). 

In the last few years, the Cosmic Microwave Background has given us an 
independent way of measuring the baryon asymmetry [3]. The relative sizes of 
the Doppler peaks of the temperature anisotropy are sensitive to ry, as illustrated 
in fig. 2. The Wilkinson Microwave Anisotropy Probe first-year data fixed the 
combination n^h? = 0.0224 ± 0.0009, corresponding to [2] 



which is more accurate than the BBN determination. Apart from the Lithium 
problem, this value is seen to be in good agreement with the BBN value. The 
CMB-allowed range is shown as the vertical band in figure 1 . 

We can thus be confident in our knowledge of the baryon asymmetry. The 
question is, why is it not zero? This would be a quite natural value; a priori, 
one might expect the big bang — or in our more modern understanding, reheating 
following inflation — to produce equal numbers of particles and antiparticles. 




(1.5) 



ri = (6.14 ±0.25) x 10 



1-10 



(1.6) 
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baryon-to-photon ratio 77 

Fig. 1. Primordial abundances versus t], courtesy of R. Cyburt 

One might wonder whether the universe could be baryon-symmetric on very 
large scales, and separated into regions which are either dominated by baryons 
or antibaryons. However we know that even in the least dense regions of the 
universe there are hydrogen gas clouds, so one would expect to see an excess of 
gamma rays in the regions between baryon and antibaryon dominated regions, 
due to annihilations. These are not seen, indicating that such patches should be 
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multipole, 1 

Fig. 2. Dependence of the CMB Doppler peaks on rj. 



as large as the presently observable universe. There seems to be no plausible way 
of separating baryons and antibaryons from each other on such large scales. 

It is interesting to note that in a homogeneous, baryon-symmetric universe, 
there would still be a few baryons and antibaryons left since annihilations aren't 
perfectly efficient. But the freeze-out abundance is 

= 10-20 (1.7) 

Tl'-y Tl-y 

(see ref. [4], p. 159), which is far too small for the BBN or CMB. 

In the early days of big bang cosmology, the baryon asymmetry was consid- 
ered to be an initial condition, but in the context of inflation this idea is no longer 
tenable. Any baryon asymmetry existing before inflation would be diluted to a 
negligible value during inflation, due to the production of entropy during reheat- 
ing. 

It is impressive that A. Sakharov reaUzed the need for dynamically creating 
the baryon asymmetry in 1967 [5], more than a decade before inflation was in- 
vented. The idea was not initially taken seriously; in fact it was not referenced 
again, with respect to the idea of baryogenesis, until 1979 [6]. Now it has 1040 
citations (encouragement to those of us who are still waiting for our most interest- 
ing papers to be noticed!). It was only with the advent of grand unified theories. 
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Fig. 3. Number of publications on baryogenesis as a function of time. 

which contained the necessary ingredients for baryogenesis, that interest in the 
subject started to grow dramatically. I have documented the rate of activity in the 
field from 1967 until now in figure 3. 



2. Sakharov's Conditions for Baryogenesis 

It is traditional to start any discussion of baryogenesis with the list of three nec- 
essary ingredients needed to create a baryon asymmetry where none previously 
existed: 

1 . B violation 

2. Loss of thermal equilibrium 

3. C, CP violation 

Although these principles have come to be attributed to Sakharov, he did not 
enunciate them as clearly in his three-page paper as one might have been led to 
think, especially the second point. (Sakharov describes a scenario where a uni- 
verse which was initially contracting and with equal and opposite baryon asym- 
metry to that existing today goes through a bounce at the singularity and reverses 
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the magnitude of its baryon asymmetry.) It is easy to see why these conditions 
are necessary. The need for B (baryon) violation is obvious. Let's consider some 
examples of B violation. 

2.1. B violation 

In the standard model, B is violated by the triangle anomaly, which spoils con- 
servation of the left-handed baryon + lepton current, 

d,J'Br.+L. = ^e^0,sW:^Wf (2.1) 

where W"^ is the SU(2) field strength. As we will discuss in more detail in 
section 4, this leads to the nonperturbative sphaleron process pictured in fig. 4. 
It involves 9 left-handed (SU(2) doublet) quarks, 3 from each generation, and 3 
left-handed leptons, one from each generation. It violates B and L by 3 units 
each, 

AB ^ AL^ ±3 (2.2) 




Fig. 4. The sphaleron. 



In grand unified theories, like SU(5), there are heavy gauge bosons and 
heavy Higgs bosons Y with couplings to quarks and leptons of the form 

Xqq, Xql (2.3) 

and similarly for Y . The simultaneous existence of these two interactions imply 
that there is no consistent assignment of baryon number to X^^. Hence B is 
violated. 
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In supersymmetric models, one might choose not to impose R-parity. Then 
one can introduce dimension-4 B-violating operators, like 

iV^R^R^abc (2.4) 

which couple the right-handed top squark to right-handed quarks. One usually 
imposes R-parity to forbid proton decay, but the interaction (12.4b could be toler- 
ated if for some reason lepton number was still conserved, since then the decay 
channels p ^ tt^v and p 7r°e+ would be blocked. 



2.2. Loss of thermal equilibrium 

The second condition of Sakharov (as 1 am numbering them), loss of thermal 
equilibrium, is also easy to understand. Consider a hypothetical process 

X -.Y + B (2.5) 

where X represents some initial state with vanishing baryon number, F is a fi- 
nal state, also with i? = 0, and B represents the excess baryons produced. If 
this process is in thermal equilibrium, then by definition the rate for the inverse 
process, Y + B ^ X,\s equal to the rate for (I2.5l i: 

T{Y + B ^ X)=T{X ^Y + B) (2.6) 

No net baryon asymmetry can be produced since the inverse process destroys B 
as fast as ( 12.5b creates it. 

The classic example is out-of-equilibrium decays, where X is a heavy particle, 
such that Mx > T at the time of decay, r = l/F. In this case, the energy of 
the final state F + B is of order T, and there is no phase space for the inverse 
decay: Y + B does not have enough energy to create a heavy X boson. The rate 
for r + B ^ X is Boltzmann-suppressed, T{Y + B ^ X) e"^^^/^. 



2.3. C, CP violation 

The most subtle of Sakharov's requirements is C and CP violation. Consider 
again some process X Y + B, and suppose that C (charge conjugation) is a 
symmetry. Then the rate for the C-conjugate process, X ^Y + B,\s, the same: 

T{X ^Y + B)=T{X ^Y + B) (2.7) 

The net rate of baryon production goes like the difference of these rates, 

JO _ _ _ 

— oc T{X -^Y + B)- T{X ->Y + B) (2.8) 



and so vanishes in the case where C is a symmetry. 
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However, even if C is violated, this is not enough. We also need CP violation. 
To see this, consider an example where X decays into two left-handed or two 
right-handed quarks, 

X ^ qlQl, X ^ qnqn (2.9) 

Under CP, 

CP: QL^ qn (2.10) 

where is the antiparticle of qn, which is left-handed (in my notation, L and R 
are really keeping track of whether the fermion is an SU(2) doublet or singlet), 
whereas under C, 

C: qL^QL (2.11) 
Therefore, even though C violation implies that 

r(X ^ qLqL) ^ r(X ^ qLqL) (2.12) 
CP conservation would imply 

T{X ^ qLqL) = r{X ^ qmn) (2.13) 
and also 

T{X ^qRqR)=T{X ^qLqL) (2.14) 
Then we would have 

r{X ^ qLqL) + r{X ^ ququ) = T{X ^ qnm) + T{X ^ ^l^l) (2.15) 

As long as the initial state has equal numbers of X and X, we end up with no 
net asynmietry in quarks. The best we can get is an asymmetry between left- and 
right-handed quarks, but this is not a baryon asymmetry. 

Let us recall exactly how C, P and CP operate on scalar fields, spinors, and 
vector fields. For complex scalars. 



C 
P 
CP 

For fermions. 



(l)^(t>* (2.16) 

(p{t, x) i(f>{t, —x) 
(l){t,x) ±(l)*{t,-x) 



C: 




P : 




CP: 





i(72i^R, i^B. ^ -i(^2iJ^*L; V^^n^^* (2.17) 

tpR{t,-x), tl!R ^ tljL{t,-x), V^7°V'(i>-^) 

ia2tp*R{t, -x), tpR -ia2ip*L{t, -x), tp ij'^tp*{t, -x) 



Baryogenesis 



13 



For vectors, 



CP 



P 



C 



A^{t,x) {A°,-A){t,-x) 
A^{t,x) {-A°,A){t,-x) 



(2.18) 



To test your alertness, it is amusing to consider the following puzzle. Going 
back to the example of X — > qq, let us suppose that C and CP are both violated 
and that 



where now we ignore the distinction between qij and q^. An astute reader might 
object that if all X's decay into qq and all X's decay into qq, with nx = n-x 
initially, then eventually we will still have equal numbers of q and q, even though 
there was temporarily an excess. To avoid this outcome, there must exist at least 
one competing channel, X ^ Y, X ^ Y, such that 



and with the property that Y has a different baryon number than qq. One can 
then see that a baryon asymmetry will develop by the time all X's have decayed. 
Why is it then, that we do not have an additional fourth requirement, that a com- 
peting decay channel with the right properties exists? The reason is that this is 
guaranteed by the CPT theorem, combined with the requirement of B violation. 
CPT assures us that the total rates of decay for X and X are equal, which in the 
present example means 



The stipulation that B is violated tells us that Y has different baryon number than 
qq. Otherwise we could consistently assign the baryon number 2/3 to X and 
there would be no B violation. 

2.4. Examples of CP violation 

We now consider the conditions under which a theory has CP violation. Gen- 
erally CP violation exists if there are complex phases in coupUngs in the La- 
grangian which can't be removed by field redefinitions. Let's start with scalar 
theories. For example 



r(X qq) ^ T{X ^ qq) 



(2.19) 



T{X -^Y)^ T{X Y) 



(2.20) 



V{X ^ qq) + T{X -^Y) = T{X ^ qq) + T{X ^ Y) 



(2.21) 



C = |5(/.|2-to2|^|2_^|^|4 

- {ii^-'+g^'' + h.c.) 



(2.22) 
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where the top Une is CP conserving, and the bottom line is (potentially) CP- 
violating. This can be seen by applying the CP transformation ^ — > </>* and 
observing that the bottom line is not invariant: 

Ccpv ^ - (/x'</'*^ + 9<l>*^ + li-c.) (2.23) 

Of course, this lack of invariance only occurs if jj? and g are not both real. 
Parametrize 

M' = |M|2e^*-, 9=\9W^' (2.24) 

We can perform a field refefinition to get rid of one of these phases, e.g. (f)^, by 
letting 

(j, ^ 6-'"^"/^ (2.25) 
so that 

Ccpv = + |5|e'(^«-'^''V^ + h.c.) (2.26) 

Now under CP, only the </)^ term changes. We can see that it is only the combi- 
nation 



</'inv = 4>g — '^(I'fj. = "invariant phase" (2.27) 

which violates CP, and which is independent of field redefinitions. If (pinv = 0, 
the theory is CP-conserving. We can also write 




(2.28) 



No physical result can depend on field redefinitions, since these are arbitrary. 
Therefore any physical effect of CP violation must manifest itself through the 
combination (pi^y. This can sometimes provide a valuable check on calculations. 

An equivalent way of thinking about the freedom to do field redefinitions is the 
following. One can also define the CP transformation with a phase (f) 
As long as there exists any value of a for which this transformation is a symmetry, 
we can say that this is the real CP transformation, and CP is conserved. 

Next let's consider an example with fermions. We can put a complex phase 
into the mass of a fermion (here we consider Dirac), 



£ = ip (id — in{cos6 + i sin Oj^ j) tp 
= ^ (id - me'^'"') ip 



(2.29) 
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[In Weyl components, the mass term has the form m{'il)Le-^^'4>R + '4'r^~^^'4'l)-] 
Under CP, 

V'^ ^ i/'^7°7^i = '07^« (2.30) 

Also, by transposing the fields, 'tpe^^^^tp = —ijj'^e^^^^j^ip*. Putting these results 
together, we find that the complex mass term transforms under CP as 

= -Vi72e*^T^7V 

= Vie-'^^^V (2.31) 

So ostensibly 6 is CP-violating. However, we can remove this phase using the 
chiral field redefinition 

^ ^ e-*(^/2)^^V (2.32) 

To get an unremovable phase, we need to consider a more complicated theory. 
For example, with two Dirac fermions we can have 

2 

C = J2^iii^- mie'^'''')^i + M(V^ie'"T^^2 + h.c.) (2.33) 

i=l 

After field redefinitions we find that 

2 

^'(^^ - ^*)^' + M(V^ie'^"-^^/^-^=/^^^^V'2 + h.c.) (2.34) 

i=l 

so the invariant phase is a — ^(6*1 + 02)- 

It is also useful to think about this example in terms of 2-component Weyl 
spinors, where the mass term takes the form 

2 

^ 4,,e^^- + /^4,ie-^fi,2 + h.c. (2.35) 

One might have thought that there are 4 field redefinitions, corresponding to in- 
dependent rephasings of ipL^i and t/jii,i, so there would be enough freedom to re- 
move all the phases. However, it is only the axial vector transformations, where 
L and R components are rotated by equal and opposite phases, that can be used 
to remove 0i. Half of the available field redefinitions are useless as far as re- 
moving phases from the mass terms. We might more accurately state the general 
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rule as "the number of invariant phases equals the total number of phases in the 
Lagrangian parameters, minus the number of relevant field redefinitions which 
can be done." 

A third example of CP violation is the famous 6 term in QCD, 

OqcdF^^F'^^^'' = l0QCDe'''"'^'F^^F^^ = -2eQCDE^ ■ B'^ (2.36) 

This combination is P-odd and C-even, thus odd under CP. Its smallness is the 
strong CP problem of QCD: if ^qcd ^ 10"^", the electric dipole moment of the 
neutron would exceed experimental hmits. The EDM operator is 

EDM = —^dnnan„^5F'^''n = dn (n^L^ ■ Eur + n^^a ■ Etil 

+i{n\(T ■ BnR — n^j^a ■ Bul)^ (2.37) 

which is mostly a ■ E, since there is a cancellation between left and right compo- 
nents in the bottom line. Contrast this with the magnetic dipole moment operator, 

MDM = ^/i„n<T^^F''''n = ju„ {n\a ■ Bur + nJjCT • Bul 

+i{'n}^& ■ Eur - n'f^a ■ Eul)^ (2.38) 

which is mostly ct .B. To verify the CP properties of these operators, notice that 
under CP, ct — > a, L <-> i?, ^ ^ —E, B ^ B. EDM's will play an important 
role in constraining models of baryogenesis, since the phases needed for one can 
also appear in the other. The current hmit on the neutron EDM is [7] 

|d„| < 3 X (2.39) 

In 1979 it was estimated that [8] 

d„ = 5 X IQ-'^^Oqcd ecm (2.40) 

leading to the stringent bound on Oqcd- 

Other particle EDM's provide competitive constraints on CP violation. The 
electron EDM is constrained [9], 

141 < 1-6 X 10"^^ e- cm (2.41) 

as well those of Thallium and Mercury atoms [10], 

\dTi\ < 9x10-25 ecm 

\dHg\ < 2x10-28 ecm (2.42) 
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2.5. More history 

I conclude this section on "Sakharov's Laws" with a bit more about the early 
history of baryogenesis. The earliest papers on baryogenesis from the era when 
GUT's were invented were not at first aware of Sakharov's contribution. Yoshimura 
was the first to attempt to get baryogenesis from the SU(5) theory, in 1978 [11]. 
However, this initial attempt was flawed because it used only scatterings, not 
decays, and therefore did not incorporate the loss of thermal equilibrium. This 
defect was pointed out by Barr [6], following work of Toussaint et al. [12], and 
corrected by Dimopoulos and Susskind [13], who were the first to use the out-of- 
equilibrium decay mechanism. None of these authors knew about the Sakharov 
paper Weinberg [14] and Yoshimura [15] quickly followed Dimopoulos and 
Susskind with more quantitative calculations of the baryon asymmetry in GUT 
models. 

The first paper to attribute to Sakharov the need for going out of equilibrium, 
by Ignatiev et al, [16] was not in the context of GUTS. They constructed a 
lower-energy model of baryogenesis, having all the necessary ingredients, by 
a modification of the electroweak theory. Their idea was not as popular as GUT 
baryogenesis, but it did make people aware of Sakharov's early paper. 

3. Example: GUT baryogenesis 

I will now illustrate the use of Sakharov's required ingredients in the SU(5) GUT 
theory [17]. This is no longer an attractive theory for baryogenesis because it 
requires a higher reheating temperature from inflation than is desired, from the 
point of view of avoiding unwanted relics, like monopoles and heavy graviti- 
nos. Nevertheless, it beautifully illustrates the principles, and it is also useful for 
understanding leptogenesis, since the two approaches are mathematically quite 
similar. 

In SU(5) there exist gauge bosons X^^ whose SU(3), SU(2) and U(l)y quan- 
tum numbers are (3, 2, — |), as well as Higgs bosons Y with quantum numbers 
(3, 1, i), and whose couplings are similar to those of the vectors. The couplings 
to quarks and leptons are shown in figure |5] We see that the requirement of B 
violation is satisfied in this theory, and the CP-violating decays of any of these 
scalars can lead to a baryon asymmetry. To get CP violation, we need the matrix 
Yukawa couplings hij,yij to be complex, where i,j are generation indices. It 
will become apparent what the invariant phases are. 

Let's consider the requirement for or Y to decay out of equihbrium. The 
decay rate is 



To — amN/^ 



(3.1) 
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Fig. 5. GUT couplings 



where a = g"^ /{Att) or ^ (y + hY /{Air), for X^^ or Y, resepctively, m is the 
mass, N is the number of decay channels, and the Lorentz gamma factor can be 
roughly estimated at high temperature as 

M ^ (1 = QT^/rr^f/^ (3.2) 
m 

considering that {E) = 3T for highly relativistic particles. The 7 factor will 
actually not be necessary for us, because we are interested in cases where the 
particle decays at low temperatures compared to its mass, so that the decays will 
occur out of equilibrium. Thus we can set 7 = 1. The age of the universe is 
r = 1 /H = 1 /F/j at the time of the decay, so we set 

Td = -H" = Vd^irr- = otmN 



— > {aNmMp/^^y'^ < m (3.3) 

where is the number of relativistic degrees of freedom at the given tempera- 
ture. Therefore we need a <C {m/Mp){y^/N). Let us first consider the decays 
of the X gauge boson. Since it couples to everything, the number of decay chan- 
nels N is of the same order as 5*, and we get 

f « 17^ (3.4) 

But unification occurs at the scale m = 10^® GeV, for values of such that 
a= 1/25, so this condition cannot be fulfilled. On the other hand, for the Higgs 
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^2 y2 ^ 
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(3.5) 



which can be satisfied by taking small Yukawa couplings. Therefore the Higgs 
bosons are the promising candidates for decaying out of equilibrium. 

For clarity, I will now focus on a particular decay channel, Y unen, whose 
interaction Lagrangian is 



UijY'^URa.ie'ji.j + h.c. 



(3.6) 



At tree level, we cannot generate any difference between the squared matrix ele- 
ments, and we find that 



IMy^e^u^l = \M 



(3.7) 



The complex phases are irrelevant at tree level. To get a CP-violating effect, we 
need interference between the tree amplitude and loop corrections, as shown in 
figure |6l 



Y 





Fig. 6. Tree plus 1-loop contributions to y — > cruji 

The one-loop diagrams develop imaginary parts which interfere with the phase 
of the tree diagram in a CP-violating manner. We can write the amplitude as 

-iM = -i {y,, + v.jFx {Ml/p") + {h*H^y),^FY {Ml/p^)) (3.8) 

where p is the 4-momentum of the decaying Y boson, 

d^g 1 1 1 



Fy = 



(27r)4 ^2 _ 4+y + te 4 



Ry + ily 



(3.9) 



and Fx is similar, but with My — > Mx and 7^ factors in between the propa- 
gators. (It will turn out that the X exchange diagram does not contribute to the 
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CP- violating interference.) In the second line we have indicated that the loop in- 
tegral has real and imaginary parts. The latter arise as a consequence of unitarity, 

-i [M{pi ^ Pf) - M*{pf ^ Pi)] = J2<>^aM{pi ^ Pa)M*{pf ^ p„)(3.10) 

a 

where the sum is over all possible on-shell intermediate states and dUa is the 
Lorentz-invariant phase space measure. lfM.*{pf Pi) = M.*{pi Pf), then 
the left-hand-side is the imaginary part of A4. The right-hand-side implies that 
an imaginary part develops if intermediate states can go on shell: 

-H <)<X) 

This will be the case as long as the decaying particle is heavier than the interme- 
diate states. 

To see how a rate asymmetry comes about, schematically we can write 

M{Y^eu) = y + yFx+yFy 

M{Y^eu) = y*+y*Fx+y*FY (3.11) 

where y = h*h'^y. Then the difference between the probability of decay of Y 
and Y goes like 

\MY^euf-\Mr^,^f = [y*{l + F;,)+y*Fp][y{l + Fx) + yFY] 

- [y{l + F*x) + yF^\ + Fx) + TFy] 
= y*yFY + y*yFY - yy*FY - yy*FY 
= [yy* - y*y]'2iiY 

= -4:lm[yy*]lY 

= -4:Imtr{h*h^yy^)Iy (3.12) 

(You can check that the Fx term does not contribute to this difference, since g 
is real.) Now, unfortunately, it is easy to see using the properties of the trace 
that ti{h*H'^yy'^) is purely real, so this vanishes. However, there is a simple 
generalization which allows us to get a nonvanishing result. Suppose there are 
two or more Y bosons with different Yukawa matrices; then we get 

IMY.^euf - \My,^su\^ = -4Y,lrntTih*BhlyBy],)Iy (3.13) 

which no longer vanishes. The function ly can be found in [17], who estimate it 
to be of order 10~^ — 10"'^ for typical values of the parameters. 



Im 
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Let us now proceed to estimate the baryon asymmetry. Define the rate asym- 
metry 

which is the fraction of decays that produce a baryon excess. The sum in the 
denominator is over all possible final states. The resulting baryon density is 

riB - Tig =^nYArABq (3.15) 

A 

where Bq is the amount of baryon number produced in each decay (Bq = 1/3 
for the channel we have been considering). Of course, we should also add the 
corresponding contributions from all the competing S-violating decay channels, 
as was emphasized in the previous section. For simplicity we have shown how it 
works for one particular decay channel. Up to factors of order unity, this channel 
by itself gives a good estimate of the baryon asymmetry. 

To compute the baryon asymmetry, it is convenient to take the ratio of baryons 
to entropy first, since as we have noted in section 1, entropy density changes in 
the same way as baryon density as the universe expands. 

ns-UB _ I^tEa^ _ 90C(3) ^ 

where we used ny = {2({3)/tt'^)T^ x 3, the final factor of 3 counting the color 
states of Y. The number of degrees of freedom is 



g* 



J2 ^1 J2 ] (spin states) 

Vbosons fcimions/ 



106.75 standard model 
160.75 SU(5) 



x(- 2 if SUSY) (3.17) 



(these numbers can be found in the PDG Big Bang Cosmology review [1]). We 
can now use jlAj to convert this to rj, by multiplying by a factor of 7. Suppose 
we have SUSY SU(5); then 

Thus we only need ^ r ^ = 6 x 10^^, a quite small value, easy to arrange since 
TA ~ -^—T2'^y (3.19) 
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where e is the CP violating phase, and we noted that Iy ^ 10~^ — 10~^. Ac- 
cording to our estimate of the size of y or ft, needed to decay out of equiHbrium, 
~ 471" X (10~^ — 10~^). We can easily make -q large enough, or even too large. 



3. 1 . Washout Processes 



The foregoing was the first quantitative estimate of the baryon asymmetry, but it 
is missing an important ingredient which reduces 77: B-violating rescattering and 
inverse decay processes. We have ignored diagrams like 




The first of these, inverse decay, was assumed to be negligible by our assumption 
of out-of-equilibrium decay, but more generally we could include this process 
and quantify the degree to which it is out of equilibrium. Any of the washout 
processes shown would cause the baryon asynmietry to relax to zero if they came 
into equilibrium. 

To quantify these effects we need the Boltzmann equation for the evolution of 
baryon number. This was first carried out by Kolb and Wolfram [18], who were 
at that time postdoc and graduate student, respectively, at Caltech. The form of 
the Boltzmann equations is 



dt 



where fx is the distribution function for particle x, J cflTj is the phase space 
integral over all final and initial particles. 



nf 



6{pI - ml) 



x(27r)4,5(4) 




(3. 



'^""^ + mnuR S j dUi {fY\M{Y ^ ne)|2 - fufe\M{ue ^ Y)f) 

+ J dUi [f.JgJMiqq ^ ue)|2 - fufe\M{ue ^ qq)\^] + 



(3.21) 



and I have ignored final-state Pauli-blocking or Bose-enhancement terms, (1±/). 
We should include similar equations for any other particles which are relevant for 
the baryon asyrmnetry (Xf^, Y, quarks, leptons), and solve the coupled equations 
numerically to determine the final density of baryons. Depending on the strength 
of the scattering, the washout effect can be important. We can estimate its signif- 
icance without solving the Boltzmann equations numerically. 
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Define the abundance Yi = rii/n-y and consider a simplified model where 
X — > hb, bb (baryons or antibaryons), with 





■^66 = 


^(i + e)rx 




rx- 


■^55 — 


i(i-e)rx 




rx- 


■*bb = 


^(i-e-)rx 




rx- 


■+55 = 


^(i + e)rx 


(3.22) 



Kolb and Wolfram found that the time rate of change of the baryon abundance 
(which we could also have called ?)) is 

Yb =Txie- e) [Fx - ^q] - ^Yb [TxYx,eq + T,,^-,-,] (3.23) 

On the right-hand-side, the first term in brackets represents baryon production 
due to CP-violating decays, while the second term is due to the washout induced 
by inverse decays and scatterings. The scattering term is defined by the thermal 

average 

^bb^bb = % {{'"^bb^bb) + {v(^bb~*bb)) (3-24) 

In the treatment of Nanopoulos and Weinberg, the effect of r^;,^^^ was neglected, 
but we can estimate it using dimensional analysis as 



2 



^bb^bb ^ /.^2 , (3-25) 



(considering s-channel exchange of the X boson), again using {E) ^ 3T. Then 
if we imagine the initial production of the baryon asymmetry happens quickly 
(by the time t = 0), followed by gradual relaxation, then 

Ysit) ~ Yb(0) e'f^"""''^* (3.26) 
Using H^l/t^ 9*T^/Mp, dt ~ {Mp/g,){dT/T^), 

1 Tdtr-. —P- I r-5- ^ rr—^ (3.27) 
J g* J lOOMx 

so the maximum BAU gets diluted by a factor ~ e-('^«'^''/(ioo^^). We thus 
need < 30Mx/Mp, hence a < 0.05, using Mx/Mp - lO""*. This agrees 
with the quantitative findings of Kolb and Turner, schematically shown in figure 
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a =10 



-3 
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B 




a= 10 



(time) 



-1 



log Mx/ T 



Fig. 7. Evolution of baryon asymmetry witli witli B-violating rescattering processes included. 

|7] Notice that this value of a is consistent with the unification value, so GUT 
baryogenesis works. 

This concludes our discussion of GUT baryogenesis. As noted earlier, it is 
disfavored because it needs a very high reheat temperature after inflation, which 
could introduce monopoles or gravitinos as relics that would overdose the uni- 
verse. Another weakness, shared by many theories of baryogenesis, is a lack of 
testability. There is no way to distinguish whether baryogenesis happened via 
GUT's or some other mechanism. We now turn to the possibility of baryogenesis 
at lower temperatures, with greater potential for testability in the laboratory. 

4. B and CP violation in the Standard Model 

During the development of GUT baryogenesis, it was not known that the stan- 
dard model had a usable source of B violation. In 1976, 't Hooft showed that the 
triangle anomaly violates baryon number through a nonperturbative effect [19]. 
In ( 12. Il l we saw that the baryon current is not conserved in the presence of ex- 
ternal W boson field strengths. However this violation of B is never manifested 
in any perturbative process. It is associated with the vacuum structure of SU(N) 
gauge theories with spontaneously broken symmetry. To explain this, we need to 
introduce the concept of Chern-Simons number. 




(4.1) 



where the current Rf^ is given by 




(4.2) 
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(note: there are a number of references [20,21] in which the factor g/i is in- 
correctly given as 2(7/3; this seems to be an error that propagated without being 
checked). This current has the property that 



(4.3) 



Chern-Simons number has a topological nature which can be seen when consid- 
ering configurations which are pure gauge at some initial and final times, to and 
ti. It can be shown that 



Ncsih) - Ncsito) 



dt 



ta 



(4.4) 



an integer, which is a winding number. The gauge field is a map from the physical 
space to the manifold of the gauge group. If we consider an SU(2) subgroup of 
SU(N), and the boundary of 4D space compactified on ball, both manifolds are 
3-spheres, and the map can have nontrivial homotopy. 

We are interested in the vacuum structure of SU(2) gauge theory. Consider 
a family of static gauge field configurations with continuously varying Ncs- 
Those configuration with integer values turn out to be pure gauge everywhere, 
hence with vanishing field strength and zero energy. But to interpolate between 
two such configurations, one must pass through other configurations whose field 
strength and energy are nonvanishing. The energy versus Chern-Simons number 
has the form shown in figure [8] 




^ N 



cs 



Fig. 8. Energy of gauge field configurations as a function of Chern-Simons number. 



Each minimum is a valid perturbative vacuum state of the theory. They are 
called n-vacua. The height of the energy barrier is 



sph 



/ 



A 



47rw 
9 



9 



2Mv 



-f 



A 



(4.5) 



.9/9 9 oiw \9 

where v — 174 GeV is the Higgs field VEV, A is the Higgs quartic coupling, 
aw — 5^/47r ^ 1/30, and the function / ranges between /(O) = 1.56 and 
/(oo) ^ 2.72. 
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't Hooft discovered that tunneling occurs between n- vacua through field con- 
figurations now called instantons (Physical Review tried to suppress that name at 
first, so the term "pseudoparticle" appears in the early literature, but it soon suc- 
cumbed to the much more popular name.) The relevance of this for B violation is 
through the relation between the divergence of the left-handed baryon plus lepton 
current to the divergence of K^^\ 

d.J'B+L = Nfd.K^ (4.6) 

where iV/ = 3 is the number of families. Integrating this relation over space and 
time, the spatial divergence integrates to zero and we are left with 

Hence each instanton transition violates B and L by 3 units each — ^there is spon- 
taneous production of 9 quarks and 3 leptons, with each generation represented 
equally. 

However, the tunnehng ampUtude is proportional to 

which is so small as to never have happened during the lifetime of the universe. 
For this reason, anomalous B violation in the SM was not at first considered 
relevant for baryogenesis. But in 1985, Kuzmin, Rubakov and Shaposhnikov re- 
alized that at high temperatures, these transitions would become unsuppressed, 
due to the availability of thermal energy to hop over the barrier, instead of tunnel- 
ing through it [23]. This occurs when T > 100 GeV. The finite-T transitions are 
known as sphaleron processes, a term coined by Klinkhamer and Manton from 
the Greek, meaning "ready to fall:" it is the field configuration of Higgs and 
which sits at the top of the energy barrier between the n- vacua [20]. A thermal 
transition between n-vacua must pass through a configuration which is close to 
the sphaleron, unless T ^ -Esph- The sphaleron is a static, saddle point solution 
to the field equations. 

Evaluating the path integral for a sphaleron transition semiclassicaUy, one 
finds that the amplitude goes like 

A ~ e-'^^P'^/'^ (4.9) 

To find the actual rate of transitions, one must do a more detailed calculation 
including the fluctuation determinant around the saddle point [21]. Khlebnikov 
and Shaposhnikov obtained the rate per unit volume of sphaleron transitions 

^ = const (^) ' ' T^e---VT (4.10) 
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where the constant is a number of order unity and mwiT) is the temperature 
dependent mass of the W boson. But this formula is only valid at temperatures 
T < £^sph- In fact, at T >, mw, electroweak symmetry breaking has not yet 
occurred, and raw (T) — 0. The Higgs VEV vanishes in this symmetry-restored 
phase, and i?sph = 0. There is no longer any barrier between the n-vacua at these 
high temperatures. 

The rate of sphaleron interactions above the electroweak phase transition (EWPT) 
cannot be computed analytically; lattice computations are required. For many 
years it was believed that the parametric dependence was of the form 

I=mtr4 (4.11) 

with c ^ 1 and a^, — g'^/i'^''^) — 1/30. This is based on the idea that the 
transverse gauge bosons acquire a magnetic thermal mass of order g^T, which 
is the only relevant scale in the problem, and therefore determines ( 14.111 ) by di- 
mensional analysis. The origin of this scale can be readily understood [25] by 
the following argument. We are looking for field configurations with Ncs ^ 1; 
so using the definition of Chern-Simons number, 

d^x{g'^FAox g^£-) - 1 (4.12) 

but the energy of the configuration is determined by the temperature, 

d^xF^r^T (4.13) 



We can therefore estimate for a configuration of size R that 
R'{g'^,o.g^A^^ ^ 1 

i?'^ ^ T (4.14) 

which fixes the size of A and R, giving R ^ l/{g'^T). Hence the estimate that 
r/V ^ (g'^T)^. However it was later shown [26] that this is not quite right — 
the time scale for sphaleron transitions is actually 5^T, not g^T, giving T/V ~ 
a^r^. Careful measurements in lattice gauge theory fixed the dimensionless 
coefficient to be [27] T/V = (29 ± 6)alT'^, and more recently [28] 

^ = (25.4 ± 2.0)alT^ = (1.06 ± 0.08) x lO^^T* (4.15) 
Ironically the prefactor is such that the old a'f^T^ estimate was good using c = 1. 
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We can now determine when sphalerons were in thermal equilibrium in the 
early universe. To compute a rate we must choose a relevant volume. We can 
take the thermal volume, l/T"^, which is the average space occupied by a particle 
in the thermal bath. Then 



(4.16) 



which must be compared to the Hubble rate, H ~ .JglT"^ jM^. At very high 
temperatures, sphalerons are out of equilibrium, as illustrated in figure |9] They 
come into equilibrium when 



(4.17) 



T 



lO^'^g* ^I'^M^ - 10"^M„ - lO^^GeV 



Hence GUT baryogenesis is affected by sphalerons, after the fact. Of course 
when T falls below the EWPT temperature ~ 100 GeV, the sphaleron rate again 
falls below H. 



log(rate) 




(time) 



Fig. 9. Sphaleron rate and Hubble rate versus time. The sharp drop in the sphaleron rate occurs at the 
EWPT. 



4.1. CP violation in the SM 

Since the SM provides such a strong source of B violation, it is natural to wonder 
whether baryogenesis is possible within the SM. We must investigate whether the 
other two criteria of Sakharov can be fulfilled. 
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It is well known that CP violation exists in the CKM matrix of the SM, 

(4.18) 





Vus 


Vub 


Vcd 




Vcb 


Vtd 


Vts 


Vtb 



Cl 




-S1S3 


S1C2 


C1C2C3 


C1C2S3 




+ S2C3e'^ 


S1S2 


C1S2C3 


C2S2S3 






+£2^36"^ 


-C2C3e'^ 



This is the original parametrization, whereas in the Wolfenstein parametrization, 
Vub and Vtd contain the CP- violating phase. Where the phase resides in the CKM 
matrix can be changed by field redefinitions. How do we express the invariant 
phase? There is no unique answer: it depends upon which physical process is 
manifesting the CP violation. However C. Jarlskog tried to address the ques- 
tion in a rather general way [29]. She showed that one possible invariant is the 
combination 

J = {mf - ml){mf ~ ml){ml - ml) 

{ml - m,l){ml - m^)(m2 - m^) K 
where ii' = 3^52^301 C2C3 sin (5 

= Im-K,,yj,y,*V}*fori^j (4.19) 

This is derived by computing the determinant of the commutator of the up- and 
down-type quark mass matrices (squared): 

J = dct [ml, ml] (4.20) 

and therefore is invariant under rotations of the fields. 

The form of J has been used to argue that CP violation within the standard 
model cannot be large enough for baryogenesis. One should find a dimensionless 
measure of the strength of CP violation, using the relevant temperature of the 
universe, which must be at least of order 100 GeV for sphalerons to be effective. 
One then finds that 



J 



10^2° (4.21) 



(100 GeV) 



12 



which is much too small to account for rj ^ 10~^". 

One might doubt whether this argument is really robust. For example, in the 
original PRL of Jarlskog [29], J was defined in terms of the linear quark mass 
matrices, which yielded a formula like (|4.19l l, but with linear rather than squared 
mass differences. In the subsequent paper it was argued that actually the sign of 
a fermion mass has no absolute physical significance, so any physical quantity 
should depend on squares of masses. However, this seems to have no bearing 
on the mathematical fact that the original linear-mass definition of J is a valid 
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invariant characterization of the CP phase. If this were the physically correct 
definition, then (14.211 1 would be revised to read J ~ lO^^*', which is in the right 
ballpark. 

A more specific criticism of the argument was given in [30], who noted that the 
argument cannot be applied to KK mixing in the neutral kaon system, coming 
from the box diagram 



s 



. s 



Farrar and Shaposhnikov point out that this CP- violating effect is not proportional 
to J, and that the relevant scale is much smaller than 100 GeV — it is the mass 
of if". The idea that J/T^^ is the correct measure of CP violation only makes 
sense if all the ratios of mass to temperature can be treated as perturbatively 
small. This is clearly not the case for the ratio of the top quark mass to the K° 
mass. Could there not also be some scale lower than 100 GeV playing a role in 
the baryogenesis mechanism? Ref. [30] attempted to construct such a mechanism 
within the SM, which was refuted by [31]. There is no theorem proving that it 
is impossible to find some other mechanism that does work, but so far there are 
no convincing demonstrations, and most practitioners of baryogenesis agree that 
CP violation in the SM is too weak; one needs new sources of CP violation and 
hence new physics beyond the SM. We will see that it is rather easy to find such 
new sources; for example the MSSM has many new phases, such as in the term 
W = IJLH1H2 in the superpotential, and in the gaugino masses. 



5. Electroweak Phase Transition and Electroweak Baryogenesis 

In the previous section we showed that B violation is present in the standard 
model, and new sources of CP violation are present in low-energy extensions of 
the SM. But the remaining requirement, going out of thermal equilibrium, is not 
so easy to achieve at low energies. Recall the condition for out-of-equilibrium 
decay, (Il4] l. 

If m ~ 100 GeV, we require extraordinarily weak couplings which are hard to 
justify theoretically, and which would also be hard to verify experimentally. 

But the EWPT can provide a departure from thermal equilibrium if it is suf- 
ficiently strongly first order. The difference between first and second order is 
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determined by the behavior of the Higgs potential at finite temperature, as shown 
in figure [TO] In a first order transition, the potential develops a bump which sep- 
arates the symmetric and broken phases, while in a second order transition or a 
smooth cross-over there is no bump, merely a change in sign of the curvature of 
the potential at = 0. The critical temperature Tc is defined to be the tem- 
perature at which the two minima are degenerate in the first order case, or the 
temperature at which V"{Q) = in the second order case. 



Fig. 10. Schematic illustration of Higgs potential evolution with temperature for first (left) and second 
(right) order phase transition. 

A first order transition proceeds by bubble nucleation (fig. fTTT l. where inside 
the bubbles the Higgs VEV and particle masses are nonzero, while they are still 
vanishing in the exterior symmetric phase. The bubbles expand to eventually 
collide and fill all of space. If the Higgs VEV v is large enough inside the bub- 
bles, sphalerons can be out of equilibrium in the interior regions, while still in 
equilibrium outside of the bubbles. A rough analogy to GUT baryogenesis is that 
sphalerons outside the bubbles correspond to B-violating Y boson decays, which 
are fast, while sphalerons inside the bubbles are like the B-violating inverse Y de- 
cays. The latter should be slow; otherwise they will relax the baryon asymmetry 
back to zero. 

In a second order EWPT, even though the sphalerons go from being in equi- 
librium to out of equilibrium, they do so in a continuous way, and uniformly 
throughout space. To see why the difference between these two situations is im- 
portant, we can sketch the basic mechanism of electroweak baryogenesis, due to 
Cohen, Kaplan and Nelson [32]. The situation is illustrated in figure [12] which 
portrays a section of a bubble wall moving to the right. Because of CP-violating 
interactions in the bubble wall, we get different amounts of quantum mechanical 
reflection of right- and left-handed quarks (or of quarks and antiquarks). This 
leads to a chiral asymmetry in the vicinity of the wall. There is an excess of 
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Fig. 1 1. Bubble nucleation during a first-order EWPT. 



Ql + 9j? relative to q^j + qj^ in front of the wall, and a compensating deficit of 
this quantity on the other side of the wall. This CP asymmetry is schematically 
shown in figure [T3] 



Sphalerons interact only with q^, not qn, and they try to relax the CP-asymmetry 
to zero. DiagramaticaUy, 




simply because there are more q^ than q^ in front of the wall. But the first in- 
teraction violates baryon number by —3 units while the second has AS = 3. 
Therefore the CP asymmetry gets converted into a baryon asymmetry in front of 
the wall (but not behind, since we presume that sphaleron interactions are essen- 
tially shut off because of the large Higgs VEV). Schematically the initial baryon 
asymmetry takes the form of figure fT4l 
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Fig. 12. CP-violating reflection and transmission of quarks at the moving bubble wall. 



CP asymmetry 



broken 
phase 



sym. phase 
wall 



Fig. 13. The CP asymmetry which develops near the bubble wall. 



If the baryon asymmetry remained in front of the wall, eventually sphalerons 
would cause it to relax to zero, because there are other processes besides sphalerons 
in the plasma which can relax the CP asymmetry, for example strong SU(3) 
sphalerons which change chirality Qs by 12 units, 2 for each flavor, as shown 
in figure [15] (See [33] for a lattice computation of the strong sphaleron rate.) 
The combination of weak and strong sphalerons would relax and B + L io 
zero if the wall was not moving. But due to the wall motion, there is a tendency 
for baryons to diffuse into the broken phase, inside the bubble. If Esph/T is large 
enough, Tsph is out of equilibrium and B violation is too slow to relax B to zero. 
This is the essence of electroweak baryogenesis. 
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Fig. 14. The B asymmetry which initially develops in front of the bubble wall. 




Fig. 15. The strong sphaleron. 



J. 1 . Strength of the phase transition 

Now that we appreciate the importance of the first order phase transition, we can 
try to compute whether it occurs or not. The basic tool for doing so is the finite- 
temperature effective potential of the Higgs field, defined by the path integral 



-/3 / d^xV^tiiH) 



- (It J d^x£[H,< 



(5.2) 



where (3 = 1/T, H is the background Higgs field, and (pi are fluctuations of all 
fields which couple to the Higgs, including H itself. Here t is imaginary time and 
the fields are given periodic (antiperiodic) boundary conditions between r = 
and /3 if they are bosons (fermions). For a compact introduction to field theory at 
finite temperature, see ref. [34]. 

To evaluate VeS, we can use perturbation theory. Since there are no external 
legs, Ves is given by a series of vacuum bubbles, 

Veff= O +(00 + ®}+... 
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The one-loop contribution has a familiar form. 



This is the free energy of a relativistic gas of bosons or fermions. Recall that the 
partition function for a free nonrelativistic gas of N particles in a volume V is 

F = -T\tlZ- 

Z = (Y.J dSe-^P'^^"^ (5.4) 

Eq. (I5.3l l is the relativistic generalization, taking into account that a fermion loop 
comes with a factor of —1. The nontrivial aspect is that we must evaluate the 
particle masses as a function of the Higgs VEV H. 

In the standard model, the important particles contributing to the thermal par- 
tition function, and their field-dependent masses, are 

top quark: mt — yH (5.5) 

( \ 



gauge bosons, W'^,W'\ B : = ^H^ 



9^ 99' 

2 



(5.6) 



V 99' 9' ) 



(One must of course find the eigenvalues of the gauge boson mass matrix.) For 
the Higgs bosons, we need to refer to the tree-level Higgs potential 



1 .V 1 



V ^ A [\H\^ - -v^j = -A + v'j (5.7) 

where the real components of the Higgs doublet are _ff = -^{(p+iXi, X2+*X3)- 
(/> is the massive component, while Xi the Goldstone bosons which get "eaten" 
by the W gauge bosons. In the effective potential, we set = so that the VEV 
is given hy H ~ ^(p, but to find the masses of the Goldstone bosons, we must 
temporarily include their field-dependence in H, take derivatives with respect to 
Xi, and then set Xi = 0- Then 

Higgs boson, H -.mfr ^ — — ^ A(30^ - v^) ^ X(6H^ - v^) (5.8) 

Goldstone bosons, Xi : m- = ^-^ = XU'^ ~ w^) = A(2i/^ - w^) (5.9) 
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We must use these field-dependent masses in the 1-loop thermal effective poten- 
tial. Unfortunately the momentum integrals cannot be done in any enlightening 
closed form, except in the high- or low-temperature limits. We will be interested 
in the high-T case, where V\-\oop can be expanded as 



_ ^ mfr2 r 2, eachrealB \ mfT r : 

Vl-loop - 2^ X I ^ ^^^^ j^.^^^ p I _ __ I ^ 



m?r2 f 2, each real B 1 m?T f 1, B 
, , , 0, F 



+ 64^ll-#-^^rU4,DiracFl+^ll^) ^^"^^^ 



where 



I -|-21n47r- = 5.408, B 
CB-21n4 5^ 2.635, F 



(5.11) 



Let's examine the effects of these terms, starting with the term of order m|T^. 
This is the source of a T-dependent squared mass for the Higgs boson, and ex- 
plains why V" (0) > at high T. Focusing on just the terms, the contribu- 
tions to the potential are 



^l-loop = H'^T'^ 



= aT^H^ + 0{H^) (5.12) 

Kree = -Xv^ + 0{H'') (5.13) 

Flitting these together, we infer that the T-dependent Higgs mass is given by 

"4(7") = -Au^ + aT^ (5.14) 

and we can find the critical temperature of the phase transition (if it was second 
order) by solving for m% = 0: 

T^ = .[^v^2^v (5.15) 
V a y 

If the transition is first order, it is due to the "bump" in the potential, which can 
only come about because of the term in Vi-ioop- The cubic term is 

Kubic = ( 3/ + l(g'+ g"f/' + OiX)) ^ -ETH' (5.16) 



where the 0(A) contributions have been neglected, since we will see that X <^ 
is necessary for having a strongly first order phase transition. 
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With the cubic term, the potential takes the form 

Vtot =mliT)H^-ETH^ + \H'' (5.17) 
and at the critical temperature it becomes 

Vtot ^ XH^ (^H ~ (5.18) 

since this is the form which has degenerate minima at iJ and H = Vc- 
Notice that the VEV at the critical temperature, Vc, is not the same as the zero- 
temperature VEV, V. By comparing ( 15.171 ) and ( 15.181) we find that 

m^T) = ^Xvl ETc = V2Xv, (5.19) 
which allows us to solve for Vc and Tc. The result is 

Tc = , (5.20) 

^aX-\E^ 

and 

E ^ 3g^ 



(5.21) 



The ratio Vc/Tc is a measure of the strength of the phase transition. It de- 
termines how strongly the sphalerons are suppressed inside the bubbles. Recall 
that 



E. 



sph 



T g \T 



(5.22) 



appears in Fsph ^ e~^°p''/^, so the bigger v/T is at the critical temperature, the 
less washout of the baryon asymmetry will be caused by sphalerons. The total 
dilution of the baryon asymmetry inside the bubbles is the factor 



(5.23) 



where tc is the time of the phase transition. Requiring this to be not too small 
gives a bound [35] 



(5.24) 
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which we can use to infer a bound on the Higgs mass [36], by combining ( I5.21l i 
and (15.241 1. Given that 

^ SM^Gf ^ 1 

the weak coupling constant is (7 — 0.65, so (I5.24l i impHes A < Sg^/Wir — 0.128, 
and therefore 

niH ^ <S2 GeV (5.26) 

This is far below the current LEP limit mn > 115 GeV, so it is impossible to 
have a strongly first order electroweak phase transition in the standard model. 

The expressions we have derived give the impression that, regardless of the 
value of the Higgs mass, the phase transition is first order, even if very weakly. 
This is not the case, as explained by P. Arnold [37]: at high temperatures, the 
perturbative expansion parameter is not just made from dimensionless couplings, 
but rather 



(5.27) 

T=0 



due to the fact that the high-T theory is effectively 3 -dimensional, and is more 
infrared-sensitive to the gauge boson masses than is the 4D (T — 0) theory. In 
the regime of large A, nonperturbative lattice methods must be used to study the 
phase transition. It was found in [38] that in the T-mn phase diagram, there is 
a line of 1st order phase transitions ending at a critical point with a 2nd order 
transition at niH = 75 GeV. For larger mn, there is no phase transition at all, but 
only a smooth crossover. This is illustrated in figure [T6l 

Most attempts to improve this situation have relied on new scalar particles 
coupling to the Higgs field. Ref. [39] considered a singlet field S with potential 

y = 2C^|i7nS'p+^2|5,|2 (5 28) 

including a large coupling to the Higgs field, so 

= /i2 + 2C^|i7p (5.29) 

If //^ is not too large {p < 90 GeV) then the cubic ltm\--^T{n'^ + 2C,^li'^fl'^ 
is sufficiently enhanced so that vjT ^ 1 for heavy Higgs. If /z^ is too large, this 
does not function as a cubic term, since it can be Taylor-expanded in . 
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75 GeV 

Fig. 16. Phase diagram for the EWPT, from ref. [38] 



H 



5.2. EWPT in the MSSM 

It is possible to implement this mechanism in the MSSM because top squarks 
(stops) couple strongly to the Higgs [40]. The MSSM is a model with two Higgs 
doublets. Hi and H2, of which one linear combination is relatively light, and 
plays the role of the SM Higgs, as illustrated in figure [TtI 




Fig. 17. Light effective Higgs direction in the MSSM (in the limit where the A" is very heavy). 



The stops come in two kinds, t l and tp., and have mass matrix 
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where TOq, mfj and At are soft SUSY-breaking parameters. If we ignore the 
off-diagonal terms (possible if we choose At tan /3 = /i) it is easier to see what is 
happening — there are two unmixed bosons coupling to H. However, we cannot 
make both of them light. This is because rnj^ gets major contributions from stops 
at 1 loop, from the diagrams shown in fig.flSl 




Fig. 18. Stop loop contributions to the Higgs 



These give 



mf ^ I m, 



. m; 



TO? 



(5.31) 



If both II and tj^ are light, it is impossible to make the Higgs heavy enough 
to satisfy the LEP constraint. Precision electroweak constraints dicatate that II 
should be the heavy one. Otherwise corrections to the p parameter, 

Ap = W 





are too large, since couples more strongly to W and Z than does 1^. 

Although we have only discussed thermal field theory at one loop, two-loop 
effects are also important in the MSSM. The t^, diagrams 

y 



y 

with gluon and Higgs exchange give a contribution of the form 




AKff = -cT^H'' In ^ 



(5.32) 



which is different from any arising at one loop, and which shift the value of v^jT^ 
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by 



Tc 2\/2A V 8A2 A 



(5.33) 



Detailed studies of the strength of the phase transition have been done using 
dimensional reduction [41], two-loop perturbation theory (see for example [42]) 
and lattice gauge theory [43] which show that negative values of (the soft 
SUSY-breaking mass parameter for iji) are required. This may seem strange; 
why would we want a cubic term of the form 

/:cubic ^ {-\ml\+y^\H^\''f^ ? (5.34) 

To understand this, we must consider thermal corrections to the squark mass as 
well, 

m}^^-\ml\+y^\H2\'' + crT'^ (5.35) 

where cr = + ^2/^(1 + sin^ /3), which can be computed similarly to the 
thermal Higgs mass: one must calculate Vcs{H ,tfi). The insertion of 1-loop 
thermal masses into the cubic term corresponds to resumming a class of diagrams 
called daisies: 



+ X5 + + + 






O 

The most prominent effect for strengthening the phase transition comes when 
— |to^| + CftT^ = 0. The negative values of to^ correspond to a light stop. 



m^^ < 172 GeV. 



It is often said that a relatively light Higgs is also needed for a strong phase 
transition, but this is just an indirect effect of eq. (15.311 1. If we were willing to 
make niQ arbitrarily heavy, we could push niH to higher values. The relation 
is [42] 



mg ^ 100 GeV exp 



(5.36) 



For example, for niH = 120 GeV, mg = 4 TeV. This is unnaturally high since 
SUSY should be broken near the 100 GeV scale to avoid fine-tuning. And it 
looks strange to have niQ ^ 
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The conclusion is that it is possible, but not natural, to get a strong enough 
phase transition in the MSSM. However it is easy in the NMSSM, which includes 
an additional singlet Higgs, with superpotential 



W = ^iHiH2 + XSH1H2 
giving a potential 



rS 



dW 


2 


dW 




+ 


dS 



soft SUSY-breaking terms 



(5.37) 



(5.38) 



Ref. [44] finds large regions of parameter space where the phase transition is 
strong enough. 

Similar studies have been done for the general 2-Higgs doublet model (2HDM) 
[45,46], with 



V 



{4hIh2 + h.c. 

1 

2' 



1 



'-MHi\ 



1 



+ HHi\''\H2\'' + \a\h\H2\'' + -X^ {H\H2f+h.c. 



-X2\H2 



(5.39) 



There are several extra scalar degrees of freedom coupling to the light Higgs in 
this model, giving a large region of parameter space where the EWPT is strong. 

A model-independent approach to increasing the strength of the EWPT has 
been presented in [47], where the effects of heavy particles interacting with the 
Higgs are parametrized by a low-energy effective potential. 



v^xwm' 



2\ ^ 



1 

A2 



(5.40) 



It is found that 1; /T > 1 is satisfied in a sizeable region of parameter space, as 
roughly depicted in figure [19] This could provide another way of understand- 
ing the enhancement of the strength of the phase transition in models like the 
NMSSM or the 2HDM, in cases where the new particles coupling to H are not 
as light as one might have thought, based on the idea of increasing the cubic 
coupling. 



6. A model of Electroweak Baryogenesis: 2HDM 

Let us now go through a detailed example of how electroweak baryogenesis could 
work. The simplest example is the two Higgs doublet model, just mentioned 
in the previous section, where it is assumed that the top quark couples only to 
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one of the Higgs fields. This assumption is useful for avoiding unwanted new 
contributions to flavor-changing neutral currents (FCNC's). 

The basic idea is as follows. There are phases in the potential ( 15.39b . in the 
complex parameters /I3 and A5. One of these can be removed by a field redefini- 
tion, leaving an invariant phase 

(6.1) 

The bubble wall will be described by a Higgs field profile sketched in figure l20l 

H,{z) = -L/i^(^)e*e.(-) (6.2) 

where 



^ sin (3 ^ 

h{z) = ^(l-tanh(l)) (6.3) 
This shape can be understood by considering the SM, where 



44 



J. M. Cline 



E 



At the critical temperature, the energy is minimized when 

dlh = 2\h{h-Vc){2h-Vc) 



(6.4) 



(6.5) 



Substituting the ansatz (16.31 1 for h into ( 16.51 ). one finds that it is a solution provided 
that 



L = 



1 1 



(6.6) 




broken 
phase 



symmetric 
phase 



Fig. 20. Shape of Higgs field in tlie bubble wall. 

We can also solve for the phases 6*;. In fact, only the difference 



(6.7) 



appears in the potential ( 15.391 ). due to U(l)y gauge invariance; it can be rewritten 



i 

= i^A,/i^ + i(A3 + A4 + A5Cos26i)/i?/i; 



(6.8) 
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where now all parameters are real, and the invariant CP phase (j) has been shown 
explicitly. (Of course, we could move it into the A5 term if we wanted to, by 

shifting e ^0-(l).)lfhi = h, then 



dV 

and c^9i + s^02 = 0, while dz{h?dzd) = (c^^ + s^^)^. These equations 
allow us to solve for 61 and 62 in the vicinity of the bubble wall. Notice that the 
trivial solution 9 = constant is prevented if ^ ^ since 

^ = ij,lhih2sin{9 + (j)) - 1-Xr,sm29hlhl (6.10) 
o9 2 

vanishes at different values of 9 depending on the value of z. In fact, as z ^ 
00, 9 —(f), while for z —00, 9 approaches a different value due to the 
nonvanishing of the quartic term in the broken phase. This guarantees that there 
will be a nontrivial profile for 9 as well as h in the bubble wall. 

The result is a complex phase for the masses of the quarks inside the wall. 
Focusing on the top quark since it couples most strongly to the Higgs, this mass 
term is given by 

-^h2iz)te'^^''H = -^h2{z)i{cose + i-f5sm92)t (6.11) 
V 2 \/2 

This means that the propagation of the quark inside the wall will have CP vio- 
lating effects. We can no longer remove 62 by a field redefinition, as we did in 
section 2, because 02 is a function of z, not just a constant. If we try to remove 
it using t e~^^'^^'^/'^t, although this removes 62 from the mass term, it also 
generates a new interaction from the kinetic term, 

iiBt iiM + ^fS6»275i (6.12) 

Any physics we derive must the be same using either description, but we will 

keep the phase in the mass term as this has a clearer interpretation. 

Now let's consider the Dirac equation, which can be written with the help of 
the chiral projection operators Pl,r = ^(1 T 75) as 



{id - thPl - m*PR)^ = 



(6.13) 
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using the complex mass 



Decompose the Dirac spinor into its chiral components as 

'Rs 



tjj = e 



-iEt 



'Xs 



(6.14) 



(6.15) 



where Xs is a 2-component spinor and s = ± 1 labels spin up and down along the 
direction of motion, which we can take to be z. The Dirac equation becomes two 
coupled equations, 



(E — isdz)Ls = mRs 
{E + isdz)Rs = m*Ls 

These can be converted to uncoupled second-order equations. 



(6.16) 



{E + isd,)-{E-isd,) 
m 



(E-isd^ — (E + isd^ 
m* 



L, = 
Rs = 



(6.17) 



We see that the LH and RH components can propagate differently in the bubble 
wall — a signal of CP violation. We can think of it as quantum mechanical scat- 
tering. 




where left-handed components incident from the symmetric phase reflect into 
right-handed components and vice-versa. One can see that as z ^ oo, m ^ 0, 
so it is consistent to have RH and LH components moving in opposite directions. 
On the broken-phase side of the wall, both components must be present because 
the mass mixes them. 
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However, quantum mechanical (QM) refelection is only a strong effect if the 
potential is sharp compared to the de Broglie wavelength. In our case, 

A^:^, L^— (6.18) 
61 mil 

while for strong scattering we need 

= ^ < 1 (6.19) 

A nih 

which is not true, except in low-energy part of the particle distribution functions. 
Thus only some low-momentum fraction of the quarks in the thermal plasma 
will experience significant QM scattering. However there is a classical effect 
which applies to the dominant momentum components with p ~ 3T: a CP- 
violating force. It can be deduced by solving the Dirac equation in the WKB 
approximation. Making the ansatz 



- A. 



,(z)e'/'p("')''"'-'^' (6.20) 



and substituting it into the equation of motion, we can solve forp(z) in an expan- 
sion in derivatives of the background fields |m(z)| and d{z). This approach was 
pioneered in ref. [48] and refined in [49]. We find that the canonical momentum 
is given by 

where 



Piz) ^po + sc 0' + 0[—] (6.21) 



PQ = signpv/l^~KM (6.22) 

is the kinetic momentum, Sc = +1 (—1) for particles (antiparticles), and the ± 
in (I6.2II 1 is + for LH and — for RH states — ^but it will be shown that this term 
has no physical significance. 

The dispersion relation corresponding to ( 16.21b is 

1 



E = y^{p±ScO'/2y^ + |m|2 - -ssc9' (6.23) 

and the group velocity is 

'dE\ po 



This is a physically meaningful result: particles with different spin or Sc are 
speeded up or slowed down due to the CP-violating effect. Notice that the effect 
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vanishes if 9' ~ or m = 0. Furthermore the effect of the ±9' /2 term in ( I6.21l i 
has dropped out, showing that indeed it is unphysical. It can be thought of as a 
being Hke a gauge transformation. One can understand this from the fact that it 
has no dependence on m. In the limit that to — > 0, the phase 9 has no meaning, 
so any physical effects of 9 must be accompanied by powers of m. The same 
can be said concerning derivatives of 9. We know that if 9 is constant, it can be 
removed by a field redefinition, hence any physical effects can only depend on 9' 
or higher derivatives. 

We can furthermore identify a CP-violating force which is responsible for the 
above effect. The Hamilton equation which gives the force is 

'"^^ - ' ss^-^ (6.25) 



dz J p 2E 2E^ 

The first term on the r.h.s. is the CP-conserving force due to the spatial variation 
of the mass, while the second term is the CP-violating force. This is the main 
result which allows us to compute the baryon asymmetry, since p appears in the 
Boltzmann equation, 

d.f P df df . 

— — I — h p — = collision terms (6.26) 

ot m oz op 

To set up the Boltzmann equations we make an ansatz for the distribution 
function of particle species i, in the rest frame of the bubble wall, 

= e/3(7(-B.+^»P.)-M0 ± 1 + ^^■^''^ 

which can be understood as follows. ^{Ei + v^Pz) is the Lorentz transformation 
of the energy when we go to the rest frame of the wall, which is moving at speed 
Vw, and is the chemical potential. The first term in (I6.27l i tells us about the 
departure from chemical equilibrium induced by the CP violating force, while 
the term 5fi is included to parametrize the departure from kinetic equilibrium. 
We are free to impose the constraint 

d^x5f^^Q (6.28) 

since /i,; already accounts for the change in the overall normalization of fi. 
The idea is to make a perturbative expansion where 



(6.29) 
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We will thus linearize the Boltzmann equations in /Uj and 6fi, and take the differ- 
ences between particles and antiparticles. The next step is to make a truncation 
of the full Boltzmann equations (BE) by taking the first two moments, 

j d^p{BE) and j d^pp^{BE) (6.30) 

The reason for taking two moments is that we have two functions to determine, 

z) and z); a more accurate treatment would involve more parameters in 
the ansatz for the distribution functions and taking higher moments of the BE. In 
any case, the truncation is necessary because the fuU BE is an integro-differential 
equation, which we do not know how to solve, whereas the moments are con- 
ventional first-order PDE's. Moreover in the rest frame of the wall, the equations 
no longer depend on time and therefore become ODE's. We can solve for 5f 
in terms of /i and further convert the coupled first-order ODE's into uncoupled 
second-order ODE's for the chemical potentials. These are known as diffusion 
equations, and they have the form 

- (AC + + Y,Vu{±^i±^i±... ) + ... = Si (6.31) 

k 

where Di is the diffusion constant, defined in terms of the total scattering rate Fj 
for particle i (in the thermal plasma) as 

Di = \^ (6.32) 

J- i 

with the thermal average (•) applied to the ^-component of the particle's squared 

velocity; Vu, is the wall velocity, usually in the range 0.1 — 0.01 depending on 
details of the model; is the rescaled chemical potential 

?i = —jT- (6-33) 
related to the particle densities by 
T ' / 1, fermions 

ni-n-, = —i\ (6.34) 
6 1^ 2, bosons 

and Tk is the fcth reaction rate for particle in which the particle i is either pro- 
duced or consumed — for example, if i + j I + m, then we would write 
r'fe(6 + Cm ~ ~ ^j); Si is the source term due to the CP violating force, 
which turns out to be 
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The source term has an asymmetric shape similar to figure |2T1 when m{z) and 
6{z) have kink-like profiles. 




Fig. 21. Typical shape of source term in diffusion equation. 



We can solve the diffusion equation using Green's functions, 



(6.36) 



where 

k± 

so 



D- 



r 6-"+^, x>o, 

k+ - fc_ \ e-'^-^, a; < 



2D 



i± Wi 



2ri:> \ 



^(z) = / dzo G(z - zo)5(zo) 



(6.37) 



(6.38) 



This tells us the asymmetry in one helicity of t quarks, say (+), which must be 
equal and opposite to that for the other helicity (— ). So far we have calculated 
the CP asymmetry near the wall. It has a shape similar to that shown in figurel22l 
The final step is to compute the baryon asymmetry. Since the asymmetry 
biases sphalerons, we get 



•Jl^ sphS sph rp2 



dt 



(6.39) 



where SFsphC is the term which causes the initial conversion of the CP asym- 
metry into a baryon asymmetry, and —cTs-phnB/T^ is the washout term, which 
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^ = |i / T = CP asymmetry 




Fig. 22. Typical sliape of tlie CP asymmetry from solving the diffusion equations. 



would eventually relax B to zero if sphalerons did not go out of equilibrium in- 
side the bubbles. In this equation, we should think of Fsph as being a function 
of t and z, which at any given position z abruptly decreases at the time t when 
the wall passes z, so the position of interest goes from being in the symmetric 
phase to being in the broken phase. (The constant c is approximately 10.) In the 
approximation that the sphaleron rate is zero in the broken phase, and ignoring 
the washout term, we should integrate ( |6.39t at a given point in space until the 
moment when the wall passes. Trading the time integral for an integral over 
z, we get 



riB 



dt- 



duB 
~dt 



dz ■ 



1 drif 



3r, 



ph 



dt 



(6.40) 



This can be generalized to the case where the washout term is not ignored, giving 
an additional factor e"^^'^/''™)'^'^^?''/^ ^ in the integrand. The entire analysis for 
the 2HDM has been carried out in [46]. 

The above derivation of the transport equations gives a clear intuitive picture 
of the physics, but more rigorous treatments have been given in [50]. 



6.1. EDM constraints 



An interesting experimental prediction is that CP violation in EWBG should also 
lead to an observable EDM for the neutron and perhaps other particles. One needs 
the invariant CP phase (f) to be in the range 10^^ — 1 for sufficient baryogene- 
sis, which leads to mixing between CP-even and odd Higgs bosons (scalars and 
pseudoscalars). Weinberg [51] notes that if any of the propagators {H2 H^*), 
{H+H+*), {H+H+*), {H+H+*), {H+H+*) have imaginary parts then CP is 
violated, giving rise to quark EDM's at one loop, through the diagram 
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t,c 



t,c 



-d 



y 

The CP-violating propagators occur because the invariant phase <j) appears in the 
mass matrix of the Higgs bosons. However this diagram with charged scalars 
requires 3 or more Higgs doublets since one of the charged fields is a Goldstone 
boson, and can be set to zero by going to unitary gauge. A similar diagram with 
neutral Higgs exchange could work, but it is suppressed for the relevant quarks 
(u and d, which are valence quarks of the neutron) because the Yukawa couplings 
are small. Moreover we should avoid coupling Hi and H2 to the same flavor of 
quarks because this tends to give large FCNC contributions. 




Fig. 23. The Barr-Zee [53] contribution to the neutron EDM. 



It was realized that two-loop diagrams give the largest contribution (see [52], 
following work of [53]), which look like figure|23] where the cross on the Higgs 
line is the CP-violating mass insertion. These contribute to the c/zrowoelectric 
dipole moment of the quark, 

CEDM = -'-gaqa^^j.G^-'q (6.41) 

where G^"^ is the gluon field strength. The quark CEDM's in turn contribute to 
the neutron EDM 

d„ = (1 ± 0.5)(0.55e di + l.le d^) + 0{dZ, d}) (6.42) 

where the factor (1 ± 0.5) is due to hadronic uncertainties. There are various 
ways of estimating the quark contribution to the neutron EDM; ( 16.421 ) uses QCD 
sum rules. See [54] for a review. 
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Ref . [46] finds that the experimental Umit on rf„ is close to being saturated for 

values of (/> which give enough baryogenesis: dn = (1.2 — 2.2) x 10~^®ecm, 
which should be compared to the experimental Umit < 3 x lO^^^ecm. 

There are extra Higgs particles with masses 300 GeV whose strong cou- 
plings to the light Higgs give rise to a first order phase transition, as we have 
discussed in the previous section. These are not CP eigenstates, so they have 
interactions which could distinguish them from the Higgs sector of the MSSM. 



7. EWBGintheMSSM 

We have already discussed how a light in can give a strongly first order EWPT 
in the MSSM. But what is the mechanism for baryogenesis? We cannot use the 
same one as in the 2HDM because in the MSSM, the Higgs Lagrangian has no 
{hIH2)^ coupling at tree level. Although such a coupling is generated at one 
loop, its coefficient is too small to get enough baryon production. 

The principal mechanism which has been considered is chargino-wino reflec- 
tion at the bubble wall. The first computation of this effect using the WKB 
approach, described in the previous section, was done by ref. [55]. (For ear- 
lier approaches, based on quantum mechanical scattering or the quantum Boltz- 
mann equation, see [56] or [57], respectively.) The WKB method has been 
carefully scrutinized and verified by [44] and a number of other papers [58]. 
Charginos/winos have a 2 x 2 Dirac mass term 




(7.1) 



where /x and m2 are complex. Denoting the 2x2 Dirac mass matrix by A^, we 
diagonalize M locally in the bubble wall using a ^;-dependent similarity trans- 
formation 

SO similarly to the top quark in the 2HDM, charginos and winos experience a 
CP- violating force in the wall, where 



mld.O^ = ±9' t iH.{z)H,iz)) (7.3) 



However the resulting chiral asymmetry in h or W does not have a direct effect 
on sphalerons. Instead, we must take into account decays and scatterings such as 
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which transfer part of the h CP asymmetry into quarks. One must solve the 
coupled Boltzmann equations for all the species to find the quark CP asymmetry 
which biases sphalerons. 

Because of the indirectness of the effect, baryon generation is less efficient 
in the MSSM than in the 2HDM. The effect is largest when ~ m^, which 
happens when m2 = fi, as shown in figure |24l taken from [59]. The figure 
shows the contours of 77 in the plane of fi and 7712, with the LEP limit on the 
chargino mass superimposed, m^± > 104 GeV. This figure was prepared before 
the WMAP determination of 77, at which time the preferred value was 3 x 10^^". 
Hence the contours only go up to 5 x 10~^° in the LEP-allowed region. It is 
possible to choose other values of the MSSM parameters to slightly increase the 
baryon asymmetry to the measured value; for example taking tan /3 < 3 helps, 
and assuming a wall thickness (which we have not tried to compute carefully, but 
is only parametrized) L ^ 6/T increases ?/. The conclusion is that all relevant 
parameters, aig{rh2fJ,), Vw, L, tan /?, must be at their optimal values to get a large 
enough baryon asymmetry. The MSSM is nearly ruled out for baryogenesis. 

In fact some additional tuning is required since the large CP-violating phases 
which are needed would lead to large EDM's through the diagram 




If the squark q is light, the invariant phase (which we can refer to as the phase of 
jj.) must be small, 0^ <, 10^^, whereas we need 6*^ ^ 1. To suppress the EDM's 
it is necessary to take the squarks masses in the range nig > TeV. Roughly, the 
bound goes like [54] 



V J 



(7.4) 




Fig. 24. Contours of constant jj X 10^ from baryogenesis in the MSSM. 



In more detail, [60] finds a bound of rriq > 3 TeV from the EDM of Hg, if 6*^ is 
maximal, as shown in figure |25] This does not conflict with the need to keep ba 
light because we only need to make df^ and (the CEDM's of the valence quarks 
of the neutron) small, while the CEDM's of higher generation quarks might be 
large. 

Assuming heavy first and second generation squarks, the largest contribution 
to EDM's can come from two-loop diagrams of the Barr-Zee type (figure |23T l, 
where the top loop is replaced by a chargino loop [61]. This contribution does 
not decouple when the squarks are heavy, and it gives stringent constraints on 
EWBG in the MSSM through the electron EDM. A heavy charged Higgs is 
needed to suppress this diagram enough to allow for maximal CP violation in 
the n parameter 

However we can loosen the constraints on the MSSM by adding a singlet 
Higgs field S with the superpotential (I5.37l i and potential ( 15.381 1. where the soft 
SUSY-breaking terms include CP violating terms of the form XASH1H2 and 
kAS^. This model has much more flexibility to get a strong EWPT and large 
CP violation [44]. We no longer need 1-loop effects to generate cubic terms 
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Fig. 25. Bound on squark masses from EDM's assuming maximal ptiase in ttie fj, parameter, from 



in the Higgs potential; already at tree level there are cubic terms H1H2S and 
S^. If S participates in the EWPT by getting a VEV, these cubic terms help to 
strengthen the phase transition, and it is possible to have Vc/Tc > 1 even with 
heavy Higgses. 

The NMSSM also provides many new CP-violating phases in the Higgs po- 
tential due to complex couplings k, fi and Ak, and the model works similarly to 
the 2HDM, generating a CP asymmetry directly in the top quark. 



8. Other mechanisms; Leptogenesis 

Unfortunately there was not enough time to cover other interesting ideas for 
baryogenesis, including the very elegant idea of Affleck and Dine which makes 
use of flat directions in supersymmetric models to generate baryon number very 
efficiently [62]. 

A very popular mechanism for baryogenesis today is leptogenesis, invented by 
Fukugita and Yanagida in 1986 [63]; see [64] for a recent review. I am not able 
to do justice to the subject here, but it deserves mention. Leptogenesis is a very 
natural mechanism, which ties in with currently observed properties of neutrinos; 
hence its popularity. Unfortunately, the simplest versions of leptogenesis occur 
at untestably high energies, similar to GUT baryogenesis. Although it is possible 
to bring leptogenesis down to the TeV scale [65], it requires a corresponding 



[60]. 
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decrease in Yukawa couplings which renders the theory still untestable by direct 

laboratory probes. 

In leptogenesis, the first step is to create a lepton asymmetry, which can be 
done without violating electric charge neutrahty of the universe since the asym- 
metry can initially reside in neutrinos. We know however that sphalerons violate 
B + L, while conserving B — L. Therefore, an excess in L will bias sphalerons 
to produce baryon number, just as an excess in left-handed B did within elec- 
troweak baryogenesis. Roughly, we can expect that any initial asymmetry Li in 
L will be converted by sphalerons into a final asymmetry in B and L given by 

Bf^-^L,, Lf^^Li (8.1) 
This estimate is not bad. A more detailed analysis gives [24, 66] 



i,SM 1 



Bf = -LA )^--Li (8.2) 

^ MSSMJ 3 ' 

since sphalerons couple to B^ + Ll and we must consider the conditions of 
thermal equiUbrium between all species. 

Leptogenesis is conceptually very similar to GUT baryogenesis, but it rehes 
on the decays of heavy (Majorana) right-handed neutrinos, through the diagrams 



* H 




+ VR_i — <^ )|cVrj 



h 

The cross represents an insertion of the L-violating heavy neutrino Majorana 
mass, which makes it impossible to assign a lepton number to vr, and hence 
makes the decay diagram L-violating. The principles are the same as in GUT 
baryogenesis. Physically this is a highly motivated theory because (1) heavy 
are needed to explain the observed vl masses, and (2) vr's are predicted by 
SO(IO) GUT'S. 

The first point is well-known; it is the seesaw mechanism, based on the neu- 
trino mass terms 

y,ji?ii,iHLj + h.c. - i {MijV^^VRj + h.c.) (8.3) 
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When the Higgs gets its VEV, the mass matrix in the basis and is 



D 



mo 
M 



(8.4) 



where the Dirac mass matrix is {mo)ij = vyij (and we must take the complex 
conjugate of ( 18.31 1 to get the mass matrix for the conjugate fields v'^, vf). This 



can be partially diagonalized to the form 







'■D 





M 



(8.5) 



so the light neutrino masses are 0{m'jj/M) ~ y^v^ /M . More exactly, they are 
the eigenvalues of the matrix v'^{yM^^y'^)ij, up to phases. We can estimate the 
size of the vn mass scale: 



M 



2 2 



y' X 10^^ GeV 



(8.6) 



where we used the tau neutrino mass, — 0.05 eV, as measured through atmo- 
spheric neutrino mixing. Since y can be small, it is possible to make M smaller 
than the gravitino bound on the reheat temperature after inflation. 

It is usually assumed that there is a hierarchy such that Mi <C M2, M3, since 
this simplifies the calculations, and it seems like a natural assumption. Then vr^i 
is the last heavy neutrino to decay out of equilibrium, and it requires the lowest 
reheat temperature to be originally brought into equilibrium. 

To find the lepton asymmetry, we need the dimensionless measure of CP vio- 
lation 



IH\ 



AH 



M. 



+ anything 



J_ M^ty)^ Ml 



(8.7) 



It can be shown that ei has a maximum value of 

^ 167r 

The baryon asymmetry can be expressed as 

rj = lO^'^eiK 



3 Mim3 



(8.8) 



where k is the efficiency factor, which takes into account the washout processes 

L\ /L Lx /t 

H, \ / 

H 



inverse 
decays 



H 



,> ^ 

AL = 2 
scatterings 



H 



AL= 1 
scatterings 
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It is interesting that leptogenesis can be related to the measured neutrino 
masses. The rate of heavy neutrino decays is given by 

r. = ^ (8.9) 

where 

mi = ^ (8.10) 

lvl\ 

is an "effective" neutrino mass; it is not directly related to the neutrino mass 
eigenvalues (which requires the transpose of rriD instead of the Hermitian con- 
jugate, and also (I8.IOI 1 is a matrix element rather than an eigenvalue), though it 
should be of the same order as the light neutrino masses. It can be shown that 
mi < mi < m3, and that the condition for Mi to decay out of equilbrium is 
satisfied if 

mi<m, = ^ — = 1.1 X lO^^eV (8.11) 

The ratio mi /m* enters into the efficiency factor k, and for mi > m*. 



k'^2x 10'^ (^ """j^" j *^^-12) 

Another interesting connection to the neutrino masses comes through the AL = 
2 scattering diagram, whose rate does relate to the actual light masses, 

rAL=2 = V ml =^T-^rn'^ (8.13) 



0.01 eV 



There is a bound on the r.m.s. neutrino mass from leptogenesis, due to this 
washout process, 

m<0.3eV (8.14) 

which intriguingly is consistent with and close to the bound on the sum of the 
neutrino masses from the CMB (see for example [67]). 

It is also possible to derive a lower bound on the lightest right-handed neutrino 
mass. Ml > 2 x 10^ GeV, assuming maximum value of the efficiency factor [68]. 
This follows from expressing the maximum value of ei in the form 
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and using Am^^^ = 2.5 x lO^^^ eV^, TO3 > Aml^^. 

These results which are consistent with the known neutrino masses (and marginally 
consistent with the gravitino bound) are suggestive that leptogenesis could indeed 
be the right theory. Unfortunately, there is no way to find out for sure through 
independent laboratory probes of the particle physics, as we hope to do in elec- 
troweak baryogenesis. On the other hand, it is much easier to quantitatively pre- 
dict the baryon asymmetry in leptogenesis than in electroweak baryogenesis; the 
Boltzmann equations are much simpler One need only solve two coupled equa- 
tions, for the heavy decaying neutrino density, and the produced lepton number. 



where z = Mi/T, D stands for the rate of decay and inverse decay, S is the rate 
of A = 1 scatterings (which can also produce part of the lepton asymmetry), and 
W is the washout terms which includes AL = 2 scatterings; see [69]. Care must 
be taken when defining the s-channel AL — 2 processes which contribute to the 
washout because when the exchanged heavy neutrino goes on shell, these are the 
same as inverse decays followed by decays, which have been separately counted; 
see [70]. 
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